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Low frequency oscillations in black holes: a spectral-timing 
approach to the case of GX 339-4 
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ABSTRACT 

We analyzed RXTE/PCA and HEXTE data of the transient black hole binary GX 339-4, col- 
lected over a time span of eight years. We studied the properties and the behavior of low 
frequency quasi periodic oscillations (QPOs) as a function of the integrated broad-band vari- 
ability and the spectral parameters during four outbursts (2002, 2004, 2007, 2010). Most of 
the QPOs could be classified following the ABC classification that has been proposed before. 
Our results show that the ABC classification can be extended to include spectral dependen- 
cies and that the three QPO types have indeed intrinsically different properties. In terms of 
the relation between QPO frequency and power-law flux, type-A and -C QPOs may follow 
the same relation, whereas the type-B QPOs trace out a very different relation. Type-B QPO 
frequencies clearly correlate with the powerlaw-flux and are connected to local increases of 
the count rate. The frequency of all QPOs observed in the rising phase of the 2002, 2007 and 
2010 outburst correlate with the disk flux. Our results can be interpreted within the framework 
of recently proposed QPO models involving Lense-Thirring precession. We suggest that type- 
C and -A QPOs might be connected and could be interpreted as being the result of the same 
phenomenon observed at different stages of the outburst evolution, while a different physical 
process produces type-B QPOs. 
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1 INTRODUCTION 

Quasi-periodic oscillations (QPOs) have been discovered in many 
systems and are thought to originate in the innermost regions of the 
accretion flows around stellar-mass black holes. Low-frequency 
QPOs (LFQPOs) with frequencies ranging from a few mHz to 
~10 Hz are a common feature in almost all black hole X-ray bina- 
ries (BHB) and were already found in seve ral sources with Ginga 



and d ivided into different classes (se e e.g _ 
1 19911 for the case of GX 339-4 and lTakizawa et al.ll 19971 for the 
case of GS 1124-68). Observations performed with the Rossi X- 
ray Timing Explorer (RXTE) have led to an extraordinary progress 
in our knowledge on properties of the variability in BHBs (see 
Ivan der Klisl l2006t IRemillard & McClintockl l2006t iBellonil I20T0I 
for recent reviews) and it was only after RXTE was l aunched that 
LFQPOs were detected in most obse rved BHBs (see Ivan der Klisl 
l2004l IRemillard & McClintockl [2006). Three main types of LFQ- 
POs, dubbed types A, B, and C, originally identifi ed in the Power 
Dens i ty Spectra (PDS) of XTE J1550-564 (se e IWiinands et"al] 
1 19991 : iHoman etal]|200ll : IRemillard et"ai]|2002h . have been seen 
in several sources. 

The systematic variations in the energy spectra of transient 
BHBs can be identified in terms of a pattern described in an 



X-ray hardness-in t ensity diagram (HID) (see IHoman et alj 200 ll ; 
IHoman & Bellonil 120051 : iBelloni et al] 120051. IBellonil l2010h . In 
many black hole candidates (BHC), different states are found 
to correspond to different branches/areas of a q-shaped HID 
pattern. Four main bright states (in addition to the quiescent 
state) have been identified in these sour ces, based on their spec- 
tral and timing properties (for a review IHoman & Bellonil 120051: 
iMcClintock & Remillardl 120061 : IBellonil l2010h . In particular, the 
analysis of the fast timing variations observed in the PDS plays a 
fundamental role in the state classification (see lHoman et al.|[2005l 
lBellonil2010l) . 

Even though the general evolution and the main transitions 
become apparent in the HID, providing a general description of 
the BHB evolution, it is not enough for detailed studies. Many ob- 
served propert ies change smooth ly throughout the basic diagrams 



served properties change smoothly throughout the basic diagrams 
(HI D, see e.g. Homan et "al] |200ll : rms vs. h ardness diagram, see 
e.g.lBellonill2010l : rms-intensity diagram, see M unoz-Darias et al. 



1201 ll) . but some do not. It is the inspection of the fast-variability 
properties which indicates the presence of abrupt variations that can 
be taken as landmarks to separate different states. In proximity of 
the HIMS/SIMS transition timing properties (in particular the ap- 
pearance of different types of QPOs in the PDS) constitute the sole 
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way to distinguish between HIMS/SIMS/HSS given the absence of 
differences in the spectral shape. 

The different types of QPOs are currently identified on the 
basis of their intrinsic properties (mainly centroid frequency and 
width, but energy dependence and phase lags as well), of the un- 
derlying broad-band noise components (noise shape and total vari- 
ability level) and of the relations among these quantities. Despite 
LFQPOs being known for several decades, their origin is still not 
understood and there is no consensus about their physical nature. 
However, the study of LFQPOs provides an indirect way to ex- 
plore the accretion flow around black holes (and neutron stars). In 
particular, their association with specific spectral states and their 
phenomenology suggests that they could be a key ingredient in un- 
derstanding the physical conditions that give origin to the different 
states. 

Several models have been proposed to explain the origin 
and the evoluti on of LFQPO in X-r ay binaries. The geome- 
try described in [Esin et al.l d 19971) and iDone et al.l d2007h allow 
to explain the spectral evolution seen in BHB s and fo r ms th e 
basis for the LFQPO mo del proposed by Ingram et alj d2009l) . 
Ingra m~& Done! (2010) and Ingram & Done (201 1), which invokes 
Lense-Thirring precession. We shall refer to this model as the pre- 
cession model. For the case of t he neutron star source 4U 1728-34, 
iTitarchuk & Osherovichl d 19991) show that LFQPO frequency is as- 
sociated with radial oscillations in the boundary layer and the break 
frequency associated to the broad band noise in the PDS is deter- 
mined by the characteristic diffu sion time of the inward motion of 
the matter in the accretion flow. lTagger&PeIiatldl999l) associate 
the existence of LFQPO in X-ray binaries to an instability occur- 
ring in the inner part of disks threaded by a moderately strong ver- 
tical poloidal magnetic field. 

In this paper, we examine the QPO parameters and compare 
them with the results of a complete spectral analysis. Since the 
presence and the properties of QPOs in BHBs are related to the 
spectral characteristics of the source, our goal is to identify the 
link between the spectral and fast timing variability properties and 
to highlight possible physical differences between the three types. 
The precession model explains both the type-C QPO frequencies 
and the presence of th e associated broad-band variability (see also 
Tagge r~& Pellad [T999I for an alternative in the case of BHs) and 
most of spectral properties. For this reason, given the success of 
this model to explain many of the observed properties we will in- 
terpret our results in the context of the precession model and we 
will attempt to extend its predictions to the three types of QPOs. 
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Figure 1. Top panel: HID for the 2002, 2004, 2007 and 2010 outburst of GX 
339-4. Each point represents a single RXTE observation. Blue stars mark 
Type-C QPOs, red squares mark Type-B QPOs and green triangles mark 
Type-A QPOs. Black dots mark all the other RXTE observations of GX 
339-4 that do not show low-frequency QPOs. Bottom panel: corresponding 
Rms-hardness diagram. Empty triangles stand for upper limits. 




1.1 GX 339-4 

GX 339-4 is a Low Mass X-Ray Binary (LMXB) harboring a > 
6M(?) accreting black hole (Hynes et al. 2003; M unoz-Darias et al.l 
l2008h . Since its discovery (Mar kert et al J 19731) , the system has un- 
dergone several outbursts, becoming one of the most studied X-ray 
transients. The intense monitoring carried out by RXTE during the 
2002, 2004, 2007 and 2010 outburst has yielded detailed studies on 
the evolution of b lack hole states throughout the outburst (see e.g. 
iBelloni et alj2005h . making the source an ideal candidate for an ex- 
tensive study on the spectral-timing properties of BHBs. Here, we 
use this rich data set to study the relation between the spectral and 
variability properties of GX 339-4 across the different outbursts. 



Figure 2. RID for the 2002, 2004, 2007 and 2010 outburst of GX 339-4. 
The symbols follow the same criteria as in Fig.fT] 



2 OBSERVATIONS AND DATA ANALYSIS 

We examined 1007 RXTE public archival observations of GX 339- 
4 from 2002, 2004, 2007 and 2010 outburst and selected for our 
analysis only observations where somewhat narrow (i.e. Q 2) 
feature was identifiable on top of peaked or power-law shaped noise 
components. The selection has been done by eye, then we fitted the 
PDS (see below) and non-significant (below 3a) detections were 
excluded from the subsequent analysis. A total of 1 15 observations 
and 117 oscillations have been considered. 

The RXTE data were obtained in several simultaneous modes. 
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STANDARD 2 and STANDARD modes for the PCA and HEXTE in- 
struments respectively were used to create background and dead 
time corrected spectra. We extracted energy spectra from PCA 
and HEXTE for each observation using the standard RXTE soft- 
ware within HEASOFT V. 6.9. Only Proportional Counter Unit 2 
from the PCA was used since only this unit was on during all 
the observations. A systematic error of 0.6% was added to the 
PCA spectra to account for residual uncertainties in the instru- 
ment calibratiorQ. For 2002, 2004 and 2007, we used only data 
coming from HEXTE/Cluster B, which was correctly working in 
that period. Since HEXTE/Cluster B encountered technical prob- 
lems at the end of 2009 [__, we decided to use data coming from 
HEXTE/Cluster A to analyze observations taken subsequently. We 
followed the standard procedure described in the RXTE cook- 
bool|f| to produce source and background spectra, using data com- 
ing from HEXTE/Cluster B to produce a preliminary background 
from which we derived a background spectrum to be used together 
with HEXTE/Cluster A spectrum. 

We accumulated background corrected PCU2 rates in the 
Standard 2 channel bands A = 4 - 44 (3.3 - 20.20 keV), B = 4 - 10 
(3.3 - 6.1 keV) and C = 1 1 - 20 (6. 1 - 10.2 keV). A is the total count 
rate, w hile the hardness was defined as H = C/B (Homan & Belloni 
l2005t) . 

PCA+HEXTE spectra were fitted with XSPEC V. 1 1 in the 
energy range 4-40 ke\|f| and 20-200 keV respectively for data 
taken in 2002, 2004 and 2007 outbursts. For HEXTE spectra pro- 
duced from data collected during 2010 we considered only the 20 
- 150 keV spectral band. The reason for this was to exclude the 
harder part of the spectra that in some cases was affected by sys- 
tematic problems due to an incorrect estimation of the background. 
For the same reason we ignored the energy range 50 - 80 keV, where 
line-like residuals can be found in the spectral fitsQ. 

We calculated unabsorbed fluxes for different spectral com- 
ponents from the best fit to the energy spectra (see Sec. 12. 2\ . We 
measured the total flux between 2.0 and 20.0 keV and the disk flux 
between 2 and 20 keV. We also measured the power-law flux in 
the 6-20 keV energy band. Since the spectral deconvolution can 
be problematic, we choose this energy interval to avoid contamina- 
tion due to the confusion between disk and power-law component 
at lower energies. All the fluxes were normalized by the Crab flux 
in the respective energy bands in order to correct the fluctuations 
due to the variations in the instrument properties. For each outburst 
we used a Crab spectrum coming from an observation as close as 
possible to the central part of the outburst. The fluxes used for the 
correction are summerized in TableQ] 

For our timing analysis, we used GoodXenon, Event and 
SINGLE Bit data modes. We used custom software under IDL and 
for each observation we produced power density spectra (PDS) 
from stretches 16 seconds long in the channel band 0-35 (2-15 
keV). We averaged the PDS and subtracted the contribution due 
to Poissonian noise ( see IZhangetalJll995r) . The PDS were nor- 
malized according to lLeahvetalJ l ll983h and converted to square 



http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmr7pcarmf- 
1 1.7 for a detailed discussion on the PCA calibration issues. 

2 http://heasarc.gsfc.nasa.gov/docs/xte/xhp_new.html 

3 http://heasarc.gsfc.nasa.gov/docs/xte/recipes/hexte.html 

4 In principle it is possible to fit RXTE/PCA spectra starting from 3 keV. 
However, especially with sources affected by low interstellar absorption, 
when fitting energy spectra residuals that cannot be explained can appear. 
Therefore we decided to fit spectra starting from 4 keV. 

5 See http://heasarc.gsfc.nasa.gov/docs/xte/whatsnew/ for details. 



Observation ID 


2 - 20 keV flux 
erg/s/cm 2 


6 - 20 keV flux 
erg/s/cm 2 


Outburst 
Outburst 


70018-01-03-00 


3.51E-08 


1.72E-08 


2002 


90129-02-01-00 


3.44E-08 


1.69E-08 


2004 


92802-03-06-01 


3.40E-08 


1.65E-08 


2007/2010 



Table 1. Crab fluxes used for the flux correction described in Sec. ff] 
Columns axe: Crab observation ID, flux measured in the 2-20 keV band, 
flux measured in the 6-20 keV band, outburst the correction was applied 
to. 
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Figure 3. Examples of type A, B and C QPOs from our GX 339-4 observa- 
tions. The centroid peak is indicated. Upper panel: Obs. 92085-01-02-06. 
Middle panel: Obs. 95409-01-15-06. Bottom panel: Obs. 70109-04-01-01). 
The Poisson noise was not subtracted. 



fractional rms ( iBelloni&HasingeJl 990 ) . The integrated fractional 
rms was calculated over the 0. 1 - 64 Hz band. PDS fitting was car- 
ried out with the standard XSPEC fitting package by using a one- 
to- one energy-frequenc y conversion and a unit response. Follow- 
ing lBellonietal] (l2002). we fitted the noise components with three 
broad Lorentzian shapes, one zero-centered and other two centered 
at a few Hz. The QPOs were fitted with one Lorentzian each, only 
occasionally needing the addition of a Gaussian component to bet- 
ter approximate the shape of the narrow peaks and to reach val- 
ues of reduced x 2 close to 1. We examined PDS in the form of 
dynamical power density spectra (DPDS), computing Fast Fourier 
Transforms of windows of data 16s long. In some cases we used 
shorter time windows to better follow the evolution of a narrow 
feature. Where transitions between different power spectral shapes 
were seen, we separated different time intervals in order to obtain 
average power spectra for each shape. 
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2.1 The QPO classification 

We have classified the QPOs following Icaselfa et al] J2004) . The 
properties that allow one to classify the QPOs are the quality factor 
(Q = fcentroid / ' FW H M) and the shape of the noise associated 
with the oscillation. Frequency does not allow a classification, as 
the characteristic frequency intervals whe re the three types of QPOs 
appear largely overl a p, nor does the rms. Case lla et al and 
Muno z-Darias et al] 1201 lh quantify the noise level associated to 
the QPOs in terms of the totajf] fractional rms (0.1-64). We could 
classify ~98% of the QPOs of GX 339-4 (see 1370. In Table[3]we 
summarize the QPOs classification and in Fig.[3]we report one ex- 
ample for each type of QPO. 

In Figure [T] we show the HID of GX 339-4 ( top panel) and 
the rms versus hardness diagram (bottom panel, see lBellonj l20ld) 
including all the RXTE observations collected during 2002, 2004, 
2007 and 2010 outbursts. Type-A, -B, -C QPOs are marked with 
green triangles, red squares and blue circles respectively. The HID 
is tracked counterclockwise, starting from the bottom right corner 
of the track. The upper and lower horizontal branches in the HID 
roughly correspond to the HIMS and SIMS, while the right and 
left vertical branches correspond respectively to the LHS and HSS. 
We will refer to the first part of the loop (from the right vertical 
branch to the left end of the top horizontal branch) as softening 
phase of the outburst and to the last part (from the left end of the 
bottom horizontal branch back to the right vertical branch) as to the 
hardening phase. 

In the following we summarize the results of the ABC classi- 
fication applied to our sample. 

Type-C QPO - In the early stages of all the four outbursts of 
GX339-4 (late LHS and HIMS, see Tab. [3] for details), two main 
components can be identified in the PDS: a strong flat-top noise and 
one or more QPO peaks. When more than one peak is observed, the 
peaks are harmonically related. The strongest and narrowest peak, 
which is usually the central one, is taken as the fundamental. When 
the identification of the fundamental remained difficult because of 
the presence of strong harmonic peaks, we followed the evolution 
of the PDS shape for the selection. In Tab. [4] we report parame- 
ters for the strongest peak in the PDS. The QPO is usually strong 
and narrow (Q 6) and the centroid frequency varies in the 0.2-9 
Hz range. Only in some cases, when the oscillations are weak and 
appear at very low frequencies, the Q-factor is slightly lower than 
10. The addition of a Gaussian component to the multi-Lorentzian 
model is required in a few observations in order to better approx- 
imate the peak shape. Type-C QPOs are observed also in the late 
stages of all the outbursts. We refer to type-C QPOs observed in 
the lower branch as type-C*. The PDS shows a noise component 
in the form of a broad Lorentzian and a QPO peak broader than at 
the beginning of an outburst. During the 2002 outburst, the type- 
C* QPOs frequencies span the 4-9 Hz range, while in the 2004 
only one Type-C* QPO (at ~ 3 Hz) is observed. In the 2007 out- 
burst Type-C* QPOs are seen in a slightly lower frequency inter- 
val (2-4 Hz). A second harmonic peak is sometimes present in the 
PDS. Even though the Type-C* QPO centroid frequency ranges, 
rms properties and Q-values are different from the case of Type-C 
QPOs, it is possible to demonstrate that the properties of the two 
kinds of QPOs are continuously connected when ordered for i n- 
creasing QPO frequency (see Sec.[3j see also lCasella et al.| [2005). 



6 We refer to the total fractional rms as to the rms measured on the whole 
PDS. 



Type-A QPO - Type-A QPOs are observed in 2002 and 2004 
during the SIMS, when the flux of the source is close to its maxi- 
mum. The PDS show a broad QPO (Q ^ 3) with centroid frequency 
between 7. 1 and 8. 1 Hz associated to a weak power-law noise. Nei- 
ther a subharmonic nor a second harmonic is observed. The PDS 
showing a Type-A QPOs have the lowest total fractional rms val- 
ues of the sample. 

Type-B QPO - They are observed in the SIMS and the oscilla- 
tions that appear in the PDS (Q ^ 6) are observable in the fre- 
quency ranges 0.8-6.4 Hz. All the type-B QPOs seen at low fre- 
quencies (i.e. below ~ 3 Hz) belongs to the lower branch in the 
HID, while all the QPOs at higher frequencies are observed in the 
upper branch in the HID. The noise seen in the PDS is weak and the 
QPO peak shape is often more similar to a Gaussian rather than a 
Lorentzian, therefore we had to combine both components (Gaus- 
sian + Lorentzian) to obtain better fits. A weak second harmonic 
is often present in the PDS. Sometimes the hint of a sub-harmonic 
appears. 

2.2 Spectral Analysis 

In order to obtain good fits and acceptable physical parameters, a 
model consisting of an exponent ially cut off power law spe ctrum 
reflected from neutral material (Magdziarz & Zdziarski 199^) was 
used (pexrav in XSPEC). The reflection parameter was left free 
to vary, while the inclination angle was fixed at 30 degrees (no- 
tice that the results only weakly depend on the inclination value 
assumed). A multi-color disk-blackbody (diskbb) was added to 
the model and a Gaussian emission line with centroid allowed to 
vary between 6.4 and 6.7 keV was further needed in order to ob- 
tain acceptable fits. The line width was constrained between 0.1 
and 1 .0 keV to prevent artificial broadening due to the response of 
XTE/PCA at 6.4 keV. The hydrogen column den sity (wabs), was 
frozen to 0.5 x 10 22 cm~ 2 ( lzdziarski et al]|2004l) . The addition of 
an iron edge, justified by the presence of the iron line, does not im- 
prove the fits significantly. In Table|4]we show the relevant spectral 
parameters for the best fits. 

Where Type-C QPOs are detected, the photon index is seen 
to rise from ~1,5 to ~2.8 and back to ~ 1.5 as a func tion of time, 
consi stently with what previously observed (see e.g. Mot ta et all 
2009). Following the loop in the HID the source becomes soft 
while approaching the HSS and subsequently becomes hard again 
going back to the LHS. As a consequence the photon index in- 
creases, remains almost constant (between ~2.6 and ~2.8) for a 
while and then decreases. When Type-A and B QPOs (during the 
softening phase) are seen in the PDS, the photon index is at its max- 
imum. When type-B QPOs are observed in the hardening phase, 
they are associated to lower values of the photon index. Indeed, 
the spectra from SIMS observations in the hardening phase always 
show systematically lower photon indices. For a detailed analysis 
of the transitions between soft and hard state in GX 339-4 dur- 
ing 2010 outburst, see Stiele et al. (Submitted). The photon index 
als o correlates wit h the L FQPO frequency, as was already noticed 
by IVignarca et al.1 120031) in the cases of GRS 1915+105, GRO 
1655-40, XTE J1550-564, XTE J1748-288 and 4U 1630-47. The 
same behavior was observed in HI 743-322 by McClintoc k et al] 
120091) . All type-C QPOs follow the same relation rather than sev- 
eral branches depending on the outburst. Type-C and type-B QPOs 
also overlap quite well covering the same photon index interval. 

The parameters associated to the iron line and to the reflection 
components do not show any clear correlations with the presence 
of the different types of QPOs and/or particular states. 
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As one can see from Table [4] not all the components of the 
model are present in all the observations. In all the observations 
where a Type-B or Type-A QPO is detected, a disk component 
(diskbb) is visible in the spectrum. When the Type-C QPOs are 
detected at hardness larger than 0.6 no disk component is observed. 
The disk appears at hardness 0.2 - 0.6 during the softening phase. 
During the hardening phase a disk component was needed only in 
some of the spectra from 2002. All the other spectral components 
(i.e. iron line, power-law and reflection components) are always 
necessary to obtain good fits. 

Examining the parameters related to the disk-blackbody com- 
ponent, it is clear that our constrains on the disc parameters are 
usually poor. Even when a soft component is clearly present in 
the spectra and is required in order to obtain good fits, it is of- 
ten only marginally significant. This is expected since the work- 
ing range of PCA (3-40 keV) allows to see only the high energy 
part of the disc black body component, above the Wien peak. It 
is also known that, even if the diskbb model provides a good de- 
scription of the thermal component, the deri ved spectral parame- 
ters should not be interpreted li terally (see e.g. lMerloni et al.l2000l 
iRemillard & McClintockl kooe). However, when the thermal com- 
ponent is dominant, the parameters can be taken as reliable. 

Disentangling the different spectral components could be 
problematic when using spectra from RXTE. If one assumes that 
the hard x-ray emission comes from Comptonization of the soft 
disk photons on hot electrons, it is known that a simple power-law 
(or a power-law-like component, such as pexrav) is not appropriate 
for the description of the hard Comptonization tail of the spectrum 
at lower energies (i.e. where the hard component overlaps with the 
soft emission from the disk blackbody). A simple power-law does 
not have the low-energy cutoff that is typical of a proper comp- 
tonization model (i.e. eqpair or compTT in XSPEC) and therefore 
could affect the real contribution of the disk. However, when us- 
ing RXTE data the adoption of a simple power-law (or power-law- 
like) component is justified because the energy range where PCA 
spectra can be analyzed (above 3 keV) does not cover this problem- 
atic overlapping zone (just below the Wien peak of the multicolor 
disk-blackbody). Therefore, a simple power-law or a power-law- 
like model such as pexrav is appropriated for the description of 
PCA spectra (see iMufioz-Darias et akllOllal and Stiele et al. Sub- 
mitted). Despite the poor constraints on the disk parameters, the 
energy spectra are well fitted by the model used and the measures 
of the disk-fluxes reported in this work are to be considered reli- 
able. This is supported by the fact that the disk flux correlates with 
the disk temperature. 



3 RESULTS 

3.1 Rms-frequency relation 

Once the QPOs of our sample were clas sified according to the ABC 
scheme, following ICasella et alj fe004h . we plotted the integrated 
fractional rms of each PDS versus the centroid frequency (see Fig. 
[4]l to probe the link between the main QPO property (the frequency) 
and the total variability of the source. Several groups of points, as- 
sociated to the Type-A, -B, -C QPOs, can be identified. 

• Type-C QPOs cover the frequency range 0.1 and 9 Hz and 
the rms range 10-35%. Type-C QPO frequency is clearly anti- 
correlated with total fractional rms. 

• Type-A QPOs form a group at frequencies in the range ~7 - 8 
Hz and rms of ~3%. 
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Figure 4. QPO centroid frequency vs. 0.1-64 Hz fractional rms. Each point 
corresponds to a different observation. Symbols correspond to QPO types: 
circles are Type-C QPOs, triangles type-A QPOs and squares type-B QPOs. 
The solid lines join for each outburst type-C QPOs during the outburst soft- 
ening, while dashed lines join type-C QPOs during the hardening phase. 
Different colors mark different outburst: black 2002, blue 2004, green 2007 
and red 2010. 



• Type-B QPOs are located at a slightly higher rms (~ 5-10%) 
in the 1-7 Hz range. 

As one can see from Figure|4] the softening phase (solid lines) 
for each outburst shows lower rms than the hardening (dashed 
lines). This property has alre ady been observed in othe r sources 
(see e.g. MAXI J1659-152, IMufioz-Darias et alj IbTQ) and can 
probably be understood in terms of a lower disk contribution to 
the emission. Only two outliers (i.e. the squares above and below 
the 5-10% of rms in Fig. [4jl can be identified in Fig. H (Obs. #10, 
#15). Obs. #10 shows the typical PDS shape of a type-B QPO, even 
though with higher rms, while Obs. #15 shows a much more noisy 
PDS. Both the points lay far both from the type-C and -B region 
in Fig. [4] However, since they follow a relation similar to type- 
B QPOs in a flux-frequency plot (see !3.3t . we tentatively classify 
those QPOs as Type-B. We notice that both the observations are 
taken in the decay phase of the outburst (2004 and 2007 respec- 
tively). 



3.2 Frequency-Power-law flux relation 

Since LFQPOs are known to be usually a sso ciated to the hard tail 
of the spectrum (see Churazov et al. 2001 and Sobolewska & Zycki 
120061 but also iRodriguez et alj [20041 and iRodriguez et alj I2008L 
who showed tha LFQPO spectra display a moving high energy 
cutoff), we started investigating the relations between LFQPO fre- 
quencies and the power-law fluxes. The frequency-power-law flux 
relation is shown in Fig. [5] We refer to power-law flux as the Crab 
corrected flux from the power-law component in the 6-20 keV en- 
ergy band (see Sec. |2.2| 7 L Different symbols represent the Type of 
QPOs and colors differentiate the four outbursts. For each outburst, 
a solid line connects in time all the QPOs from the first type-C to the 
first type-B QPO. The typical time interval between the last type-C 



7 The result depend only very weakly on the energy band chosen. We per- 
formed the powerlaw-flux measure also in the 2-20 keV and the getting 
almost identical results. 
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Figure 5. QPO centroid frequency vs. 6-20 keV powerlaw flux, normalized to the Crab flux. Colors indicate the outburst: black stands for 2002, blue for 2004, 
green 2007 and red 2010. Symbols correspond to QPO types: stars are Type-C QPOs, triangles type-A QPOs and squares type-B QPOs. The solid lines join for 
each outburst type-C QPOs and the first type-B QPO detected after the disappearing of the Type-C QPO. The empty circles mark the type-C QPOs observed 
immediately before the appearance of a type-B QPO. Each point in the plot represents an entire RXTE observation in which a QPO was detected, apart from 
the cases in which a switch between two different types of QPOs was observed (Obs. #l/#26, Obs. #4/#27-#28). 



and first type-B QPO is ~ 1 day, even though sometimes other types 
of PDS are observed in between, especially if the transition period 
is not densely observed. When this happens, it might be possible 
to miss the appearance of a transient type-B QPO and see a type- 
A QPO after a type-C QPO. In Figure 13 the different QPO types 
follow clear and separate relations as function of the hard flux. 

• Type-C QPOs (stars) lie on the right part of the diagram. The 
points trace well-defined tracks at different flux levels for each out- 
burst (see the solid lines in the plot). The whole 0.1 - 9 Hz fre- 
quency range is spanned. Each track roughly correspond to the late 
LHS and HIMS observed in the softening phase of each outburst. 

• Type-C* QPOs follow tracks on the left part of the diagram. 
Those QPOs were observed in all the outbursts. Differently from 
what happens in the right part of the diagram, points belonging 
to different outburst span different and smaller frequency ranges. 
However, the maximum frequency reached is consistent with the 
softening tracks. We ascribe the fact that no QPO appear at lower 
frequencies to the count rate being very low. 

• Type-A QPOs (triangles) cluster on a quite narrow frequency 
and flux range, close in frequency to the last Type-C QPOs seen 
before the transition to the SIMS, but at slightly lower fluxes (see 
Fig. EJl. Type-A QPOs always appear in time after the detection of 
a type-B QPO. No Type-A QPO is observed in the left part of the 
plot, i.e. during the hardening phase at the end of the outburst (see 
alsoFig.CQ.S 

• Type-B QPOs (squares) are sharply correlated with the 
power-law flux and the relation between frequency and powerlaw- 
flux is well described by a power-law of the form y — Ax B + C 



(where A = 19.4(8) , B = 0. 1 8(6), C= -6. 1(2), see Fig.©. This corre- 
lation holds for a large range of flux, showing that these oscillations 
frequencies depend directly on the hard X-ray flux. However, un- 
like type-C QPOs, the points do not follow a clear path as a function 
of time. No oscillations are seen in a given flux range in the middle 
of the plot. As happens for type-C and type-C* QPOs, this is due 
to the fact that also the SIMS (where type B QPOs are observed) is 
crossed two times, at either high or low fluxes. 

Since type-B QPO are found in a small hardness range, it 
might be argued that they behave like type-C QPOs when observed 
in a small hardness range. For this reason we checked whether type- 
C QPOs show a behavior similar to type-B QPOs once grouped in 
subsamples selected as a function of the hardness. We divided type- 
C QPOs in six subsamples and for each group we plotted the QPO 
centroid frequency as a function of the power-law flux. The result 
is shown in Fig. [7] 

Type-C QPOs' frequencies only show a weak anti-correlation 
with the power-law flux, especially at high frequencies. When type- 
C QPOs' range overlaps the hardness range where type-B QPOs are 
found (red points in Fig.[7]in the hardness range 0.2-0.3) the (weak) 
correlation that they follow is exactly opposite to the one shown by 
type-B QPOs. This fact further strengthens the difference between 
type-C and B QPOs. 

3.3 Frequency-Disk flux relation 

In Fig. [8] we show the relations between frequency and disk flux. 
Symbols and colors follow the same criteria of Fig. [5] We re- 
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Figure 6. Type-B QPOs frequency as the function of the power-law flux. 
The correlation is well described by a power-law of the form y = Ax B +C, 
where A = 19.4(8) , B = 0.18(6), C= -6.1(2) 
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Figure 7. Type-C QPO centroid frequency versus power-law flux. Differ- 
ent colors marks the different hardness ranges in which the QPOs where 
divided. From the blue to the red the ranges are: 0.7-0.8, 0.6-0.7, 0.5-0.6, 
0.4-0.5, 0.3-0.4, 0.2-0.3. 



fer to disk flux as the Crab corrected flux coming from the disk- 
blackbody component in the 2-20 keV energy range. Since we 
could associate a measure of the disk flux only to a subsample of 
type-C QPOs, not the all of them are present in this plot. 

Most of the points trace out a well-defined track and the fre- 
quency clearly correlate with the disk flux. Those points correspond 
to all the QPOs seen during the upper branch of the four outbursts. 
Notice that for three of the four outburst (2002, 2007 and 2010 out- 
burst) the upper branch in the HID loop is the same, while during 
the 2004 the upper branch is observed at a lower flux level. Note 
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Figure 8. QPO centroid frequency versus soft flux. The symbols follow the 
same criteria as in Fig. [5] 



that this track includes type-C as well as type-A and B QPOs. Type- 
A QPOs are located in correspondence of the highest disk fluxes, 
while type-B QPOs cover approximately the same flux range of 
type-C QPOs. 

We identify also a number of outliers, which correspond to 
two branches at different flux levels, forming other two tentative 
correlations. The clearest of the two is formed by all the type-C 
QPOs observed during the hardening phase (black points in the left 
upper corner of Fig. [8}. Only during 2002 outburst it was possible 
to measure the disk flux during the lower branch. We ascribe this to 
the disk being fainter and/or colder during other outbursts than in 
2002. 



3.4 Association of Type-B QPOs with local peaks in the light 
curve 

In Fig. [9] we plot sections of the light curves (PCU data, 2-20 keV) 
of the four outburst of GX 339-4. From this figure it is clear that 
most of the type-B QPOs are found at times of local peaks in the 
light curve. In all the cases where different types of QPOs are ob- 
served within a short time interval, they follow a precise count rate 
segregation: type-B QPOs at highest count rates, Type-C QPOs im- 
mediately below and type-A QPOs are found at lower count rates. 
A relation QPO-type/cou nt rate therefore seem s to exist, although 
it is different from what ICasella et alj (2005) observed in XTE 
J 1859+226, where type-A QPOs are seen at highest count rates and 
type-B and -C appears below, even though they are not clearly sep- 
arated in count rate. We notice that in all the outbursts, there is 
always a B at lower flux than some C in the same outburst, there- 
fore the segregation in count rate in not absolutely true all along 
the outburst, but only for certain intervals (see Fig.|9j. We also no- 
tice that outbursts 2002 and 2004 had a different initial evolution 
in comparison to 2007 and 2010, where the count-rate peak was 
reached after a monotonic rise. 

Despite a difference segregation in count rate, also in XTE 
J 1859+226 type-B QPOs are always found at hardness higher than 
that of type-A QPOs and lower than that of type-C QPOs. In the 
case of XTE J 1859+226 a certain overlap was observed, while there 
is no overlapping in the case of GX 339-4 during a single outburst. 
A similar correlatio n between QPO typ es and hardness was found 
in XTE J 1550-564 drloman et al] l l200lh ). 
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Figure 9. Light curves of the 2002, 2004, 2007, 2010 outburst of GX339-4. 
Each point represents one entire RXTE observation. The symbols mark the 
different kinds of QPOs following the criteria used in Fig. [JJ The Type-B 
QPOs appear associated with count-rate peaks. 



3.5 Timing and spectral evolution 

In four observations of our sample the PDS show rapid transitions 
between different shapes. In all cases the transitions involve type-B 
QPOs. 

In Fig. [To] and [TT] we show two examples of different behav- 
iors, for Obs. #l/#35 and #4/#36-#37 respectively. In the case of 
Obs. #l/#35 (Obs. ID 70109-01-07-00, Fig.[l0j, a type-A QPO (~7 
Hz) is present in the first part of the observation, when the observed 
count rate was low. In the second part the light curve shows a net in- 
crease in count rate and si multaneously the ons et of a type-B QPO 
(~6 Hz) is observed (see iNespoli et ai1l2003l for details). At the 
same time, an increase of the hard flux (from ~ 11% to ~13% of 
the total flux) is observed, as the variation in hardness suggests (see 
Tab. [3}. 

During Obs. #4/#36-#37 (Obs. ID 70108-03-02-00, second 
PCA orbit, Fig. II It a similar situation can be observed. In the first 
orbit the PDS shows a type-A QPO (~7 Hz). In the second orbit the 
source count rate dropped abruptly from ~2200 counts/s to ~2100 
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Figure 10. Upper panel: Light curve for Obs. #l/#35. The red line marks 
the light curve interval where a type-A QPO was detected. The black line 
marks the interval wh ere a type-B QPO w as visible. Lower panel: PDS for 
the two time intervals INespoli et alj2003l 



counts/jf| in few seconds. A type-B QPO (^5.6 Hz), that was visi- 
ble in the first part of the observation disappears leaving place to a 
type-A QPO (~7 Hz), observable until the end of the interval and 
during the complete third orbit. Analogous to the previous case, a 
variation in flux takes place. When the type-B QPO disappears, the 
power-law flux is seen to decrease abruptly (from ~ 15% to ~11% 
of the total flux in ~ls). The frequency of the type-A QPO before 
the appearance of the type-B QPO and and after its disappearance 
in the light curve is the same. 

Two other cases are found in Obs. #5 and #30 (Obs. ID 701 10- 
01-47-00 and 95409-01-19-00), where, in correspondence of a rise 
in the count rate, a type-B QPO takes the place of power-law- 
shaped noise (i.e. no type-A or C QPO is observable before the 
onset of the type-B QPO). In all the mentioned cases it is clear that 
spectral differences can be very subtle, much more than the timing 
changes. 

Similar fast transitions between different types of QPOs have 
already been observed in GX 339-4 dMivamotq_& Kitamotol 1991 



and in other sourc es, such as XTE J 185 9+226 dCasella et alj|2004l) 
and GS 1124-68 dTakizawa et ai]| 19971) . For XTE J1859+226, the 
type-B QPO seems to be associated to a flaring behavior. However, 
also in this source the type-A QPO is always seen at slightly higher 
frequencies than the type-B QPOs. 

As it happens for type-A/B QPOs, direct switching from/to 
type-C/t ype-B can be observed in GX339-4 as in other sources 
(see e.g. lMivamoto et alJl99lLrTakizawa et aljl997l) . However, the 
switch between type-B and -C QPOs is not as sharp as in the case of 
type-B and type-A QPOs: the transition between the two types usu- 
ally comes with a complex behavior and Type-B QPOs appears in 
correspondence to peaks in the light curve (occurring at timescales 
of few seconds), consistently with what described in Sec. 13.41 and 
type-C QPOs are seen where the count rate drops. For this reason, 
transitions between/from type-C and type-B QPOs are not easily 
detectable and are worthy of a more detailed analysis that is be- 
yond the scope of this work. For a more detailed study, see Homan 
et al. (in prep). 



Note that the counts are taken from EVENT mode data and they are not 



normalized per number ofPCUs on. 
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Figure 11. Upper panel: Light curve for Obs. #4/#36-#37 (second orbit). 
The red line marks the light curve interval where a type-A QPO was de- 
tected. The black line marks the interval where a type-B QPO was visible. 
Lower panel: PDS for the two intervals. 
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Table 2. Summary of type-A, -B and -C QPOs properties in GX 339-4. (*) 
The bracket values correspond to the hardening phases. 



4 DISCUSSION 

We analyzed RXTE/PCA and HEXTE data collected over eight 
years of observations of the transient BHB GX 339-4 to study the 
properties and the behavior of LFQPOs. 115 out of 117 oscilla- 
tions could be classified into the three main types (A, B, C). The 
coherent scenario we constructed can be compared to that of other 
systems (such as XTE 1859+226 and XTE J 1550-564) for which 
the ABC classification has been performed. Different properties 
and relations emerge from the analysis, allowing a better charac- 
terization of the different types of QPOs. Our results confirm that 
the ABC classification can be extended to include spectral depen- 
dencies. The main parameters of the different types of LFQPOs 
observed in GX 339-4 are summarized in Table|2] 

Type-C QPO : this oscillation is found in LHS and HIMS (see 
Bellonj| |2010l and Homan et al 2011 in prep.). It is observed in 
both the softening and hardening outburst intervals, even though 
during the hardening it results weaker. It appears at hardness 
values in the 0.2 - 0.8 range and over a large frequency range (0.2 
- 9 Hz). Its centroid frequency rises as function of the time as the 
source undergoes the softening phase and decreases as the source 
hardens during the decay phase. All type-C QPOs are spread above 
a variability level below which only t ype- B and type-A QPOs 
are o bserved (see Casell aet alj 1 20041 and iMunoz-Darias et ail 
|20T ]]). The observed 0.1-64 Hz rms has values between 10 and 
35% . Those are p ositively correlated with the hardness ratio (see 
e.g. iBellonil l2010h and negatively correlated to the frequency 



(see Fig. |4j- Type-C QPOs form a clear but complex pattern in 
a frequency versus power-law-flux plane. The frequency type-C 
QPOs correlates with the disk flux and form different branches in 
a frequency versus disk-flux plane, corresponding to the softening 
and hardening phases. Type-C QPO frequencies correlates well 
with the hardness, stressing a clear dependence on the spectral 
shape. 

Type-A QPO: this QPO is usually observed in the SIMS, which 
is indeed defined on the basis of the appearance of type-A and -B 
QPOs. It is found in a narrow hardness (0.20 - 0.22), frequency 
(6.5 - 8.0 Hz) and rms (~ 2 - 3%) range. The frequency at which 
they are found is always very close to the frequency of the last 
type-C QPOs observed before the transition to the SIMS. We 
observed type-A QPO only during the softening phase (i.e. along 
the upper horizontal branch in HID). The lack of this type of QPO 
during the hardening phase can however be ascribed to the lower 
statistics, as the feature is weak and broad. This is the QPO type 
that is found to be associated to the lowest total fractional rms 
in our sample. In a frequency versus power-law-flux plot, type-A 
QPOs appear to be grouped close to the high-frequency end of the 
tracks defined by the type-C QPOs and are found around the same 
frequency of the type-C QPOs observed close to the transition to 
the SIMS. 

Type-B QPO: its presence defines the SIMS. The frequency 
range where it is observed is 0.8 - 6.4 Hz. In addition, all type-B 
QPOs are observed at lower frequencies with respect to the type-C 
QPOs observed just before the transition to the SIMS. As for type- 
C QPOs, these oscillations are seen in both the softening and hard- 
ening outburst phase. Total fractional rms and hardness values are 
lower than in the case of type-C, but higher than for type-A QPOs, 
ranging in the intervals 5 - 10% and 0.2 - 0.3 respectively. It is no- 
ticeable that between 5 and 10% bro ad band total fractional rm s 
only type-B QPOs are observed (see Munoz-D arias et al]|201 lb . 
This type of QPO follows a sharp frequency-flux correlation for 
both disk and powerlaw flux (Fig. [5] and [8). 

From Fig.|5]one can see that type-B QPOs follow a clearly different 
path compared to the other classes, suggesting the presence of an 
intrinsic difference. 



It is interesting to compare some of the properties re- 
ported here for GX 33 9-4 with those observed in other sources. 
iMcClintock et al.l2009l performed a detailed spectral analysis and a 
supporting timing analysis on the BHC H1743-322 and compared 
their r esults with the BHC XTE 11550-564, while ISobczak et al] 
(2000) compared the properties of XTE 11550-564 with GRO 
11655-40. These three sources showed a behavior similar to GX 
339-4 in the frequenc y-disk flux plane, spanning a larger flux range. 
ISobczak et ail 12000) also investigated the relations between fre- 
quency and powerlaw flux, finding opposite relations for GRO 
11655-40 and XTE J1550-564. The frequency/poweiiaw-flux cor- 
relation in GX 339-4 is similar to that of GRO 11655-40, but the 
frequency range covered by the two sources do not overlap, mak- 
ing impossible a direct comparison. 

The relation between QPO frequency and power-law photon 
index is the same fo r GX 339-4 arid H1743-322, XTE 11550-56 4 
and GRO 11655-40 JSobczak et alj|200d : lMcClintock et al]|2009l). 
follow ing the original correlation observed by IVignarca et all 
(2003). As for the frequency- total rm s plane, which was orig- 
inally presented for XTE 11859+226 dCasella et all ( l2005h ). the 
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same relation is present also in H1743-322 and XTE J1550-564 
jMcClintock et alj2009h . 



4.1 Similarities and differences: a common origin for 
QPO-types? 

The different types of QPOs often show similar/compatible proper- 
ties (eg. centroid frequency range, QPO profile, quality factor val- 
ues), some of which suggest that there could be a common origin 
for the different classes. However, there are systematic differences 
that cannot be ignored. 

As is clear from Fig. [4] a stringent relation between the total 
(0.1 - 64 Hz) fractional rms values and QPO-type exists and the 
different types of QPOs correspond to different and well separate d 
rms ranges (see lCasella et alj[2004l and lMunoz-Darias et alj|20nh . 
At the same frequency two or sometimes even three different types 
of QPOs can be seen (not simultaneously) depending on the vari- 
ability level at which the system is observed. 

Despite the clear separation in rms, type-C and type-A QPOs 
follow a similar hard-flux/frequency relation. In addition, type-A 
QPOs and the type-C QPOs observed just before the HIMS/SIMS 
transition show very similar frequencies (see [4] and Tab. [3j, while 
type-B QPOs are systematically found at lower frequencies. Un- 
fortunately type-A QPOs constitute only a small part of our entire 
sample of QPOs (8 out of 1 17), therefore it is not possible to com- 
pletely exclude a bias in the frequency association due to the small 
number of detections. 

It is clear that Type-A and -C QPOs are significantly differ- 
ent timing features, in particular for what concerns the broad band 
noise associated to the two types of oscillations (strong broad band 
noise in the case of type-C QPOs and weak power-law or weak 
peaked noise in the case of type-A QPOs). However, this is not 
enough to rule out the possibility that type-A and type-C QPOs 
share a common physical origin. 

Here we discuss the evol ution of the LFQP Os in the frame- 
work of the model proposed bv lDone et alj J2007h . that also suggest 
possible explanation of the simultaneous spectral transition seen 
during an outburst. In the geometry that these authors assume, the 
outer accretion flow takes the form of a cool, geometrically thin, 
optically thick accretion disk truncated at some radius, which is 
larger than the last stable orbit (in the HIMS). The inner accretion 
flow, instead, forms a hot, geometrically thick, optically thin con- 
figuration. The inward movement of the truncation radius with in- 
creasing mass accretion rates within the hard states gives a physical 
basis for the hard-soft transition when the disk finally reaches the 
last stable orbit. The inner disk radius evolution gives also a possi- 
ble origin for the LFQPOs and to the associated noise observed in 
the power density spectra. In this scenario LFQPOs (and in partic- 
ular ty pe-C QPOs) arise f rom the vertical Lense-Thirring preces- 
sion dStella & Vietri|[l998h of the misaligned inner hot flow, while 
the broad band noise a rises from propagation of Magneto Rota- 
tional Instability (MRI. lBalbus & Hawlevl|l99llFI ) fluctuations of 
the same hot flow (se e llngram et alj|2009l llngram & Donell2010l . 
Ingram & Done 2011 for details j3 The QPO frequency depends 



9 Notice that MRI is not able to produce the jets observed during the hards 
states, therefore a different/additional mechanism must be also at work dur- 
ing the hard state. This is also valid during the transition to the SIMS where 
relativistic ejections can be observed. 

10 Fo r alt ernative models on the origin of LFQPOs see e.g. 



( 2001),[3 Titarchuk & Osherovich 1 1999), Tagger & Pellat 
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on the truncation radius and also on the optical depth of the in- 
ner hot flow that produces them, as well on the inclination of the 
system (see also Ho man et al.1 120051) and to the optical depth of 
the precessing inner hot flow. In this scenario, the fact that type- 
A QPOs and the type-C QPOs observed close to the transition are 
seen at the same frequency, suggest either that they are produced 
in a region at similar radii or in a medium at the same physical 
conditions (i.e. same optical depth). This second possibility is sup- 
ported by the clear difference in broad-band noise level in the PDS, 
that is thought to be related to MRI. The type-C QPO width sets 
the typical timescale for the QPO signal to remain coherent. Since 
the QPO is not always present in the light curve, there also ex- 
ists an excitation timesca le responsible for the QPO appearance 
dLachowicz & Dond feoiO). The QPO can be broadened if the ex- 
citation timescale becomes shorter than the QPO timescale, as the 
excitation presumably interrupts the coherent QPO light curve by 
phase randomizing it. In the framework set b y the global Lense- 
Thirring precession of the hot flow mo del {Fragile et alj 1 20071 : 
llngram et at]|2009l ; llngram & Donell201ll) the QPO could be trig- 
gered by turbulence, resulting in a vertical kick applied by the bend- 
ing waves propagating in the outer accretion flow to the precessing 
inner hot flow at random phase. If more than one kick occur within 
the QPO timescale, the QPO coherent light curve is interrupted by 
a random phase shift. As a result, the type-C QPO would evolve 
in a broader and fainter feature (Chris Done, private communica- 
tion). Thus, type-A and -C QPOs could be the result of the same 
physical process, i.e. Lense-Thirring precession. However, given 
the undeniable link between the type-A QPOs and very weak noise, 
it remains a fact that whatever process broadens and weakens the 
QPOs, should be also responsible for the collapse of the noise. 

4.2 The peculiar case of the type-B QPO 

Type-B QPOs show properties that differentiate them from the 
other two classes. Analyzing the frequency-hard flux relation and 
the frequency-rms relation, it is evident that there is a clear discon- 
tinuity within the pattern defined by type-B QPOs and other types 
of QPOs. While all the QPO frequencies seem to correlate with 
the soft flux, only type-B QPOs show a sharp correlations with the 
hard flux. This is remarkable because when type-B and type-A/-C 
QPOs are seen at similar hardness (type-C QPOs observed just be- 
fore the transition and all type-A QPOs), there are no differences 
in the spectral shape, as one might deduce from the HID, but only 
in flux (see below). In addition, type-B QPOs are systematically 
found at lower frequencies with respect to the last type-C QPOs 
and type-A QPOs ( see g] and [5}. 

Type-B QPOs can be transient (i.e appear/disappear in few 
seconds and are ob servable only for short periods, see also 
Takiza wa et al. 1 1 1997b and vary significantly around their cen- 
troid frequency, w ith a characteristic time scale of ~10s (see 
iNespoli et afll2003h . Also type-C QPOs can appear and disappear 
in few seconds, but they remain observable for long periods and 
can be easily followed in their frequency evolution during the hard- 
to-soft or soft-to-hard transition. Because of its intrinsic faintness, 
type-A QP Os cannot be follow ed as can be done for type-B and -C 
QPOs, (see lNespoli et alj|2003b . 

A noticeable peculiarity of type-B QPOs is the association to 
flux peaks. This is particularly evident for a direct switch from or 
to a type-B QPO (see Sec. 13 -5 1 ■ The association of type-B QPOs 
with increases in the count rate can be seen both in the total light 
curve (i.e. peaks in th e count rate observ ed in the total light curve 
of the source, see also lFender et alj2009b and on shorter timescales 
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(i.e. when sudden increases in count rate take place during a single 
RXTE pointing). Under the assumption that the count rate tracks 
the accretion rate, type-B QPOs would be related to increases of 
the local accretion rate. Alternatively, the increase in the count rate 
associated to those QPOs could be related to the presence of a jet 
component that would contribute to the hard emission. A third pos- 
sibility is that type-B QPOs might occur simultaneously to sudden 
changes in the geometry or radiative efficiency, which would pos- 
sibly cause variations in the relative contribution of the emitting 
component to the spectrum and in the flux. 

Starting from the precession model and following a reason- 
ing similar to that described above, we argue that type-B QPOs are 
either produced in a different region located at larger radii (with 
respect to the region where the type-C QPOs close to the transition 
and type-A QPOs might originate, to match the lower frequencies 
observed) or coming from a modulation operated by a medium with 
different physical properties. In the first case, it is necessary to find 
a process different from the vertical Lense-Thirring precession, that 
would be able to produce modulations at larger radii. This hypoth- 
esis is supported by the fact that when fast switches from/to type- 
A/-B QPOs are observed, the type-A QPOs is consisten t with being 
still p resent when the type-B QPOs appears (see also lNespoli et al.l 
2003). This is also valid for type-A and -B detected in separated ob- 
servations. Such a property suggests that two different, eventually 
simultaneous mechanisms, might be responsible of the production 
of type-C/A QPOs and type-B QPOs. For the second case, a pro- 
cess able to trigger fast transitions in the physical properties of the 
plasma would be needed. This hypothesis is supported, as for type- 
A QPOs, by the fact that the transition to type-B QPOs is associ- 
ated to the significant difference in the PDS broad band noise level. 
There is also a third possibility: type-B QPOs could be produced 
in the same region where type-C QPOs come from, but thanks to a 
different pecession mode. In the precession model the QPO arises 
from the surface density-weighed Lense-Thirring-precession over 
the inner hot flow. A sudden change in the surface density profile 
- for example due to a jet ejection from the very inner regions of 
the accretion flow - would result in a different weighing and also a 
different precession frequency. 



However, the characteristic properties of type-B QPOs suggest that 
they could be the effect of a physical phenomenon different from 
the Lense-Thirring precession and possibly somehow related to the 
transition/jet ejection mechanism. 
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5 CONCLUSIONS 

The large amount of RXTE observations of GX 339-4 in the past 
eight years allowed us to analyze the spectral and temporal behav- 
ior of the source over four outbursts. We considered all the obser- 
vations where a low frequency oscillation was observed and per- 
formed a complete spectral and timing analysis. Almost all the os- 
cillations observed in the PDS could be classified following the 
ABC classification and the three types of QPOs display differ- 
ent dependences on the spectral and timing parameters, further 
strengthening their intrinsic differences. 

We conclude that type-B QPOs show properties that clearly 
differentiate them from other types of QPOs. Their frequencies 
clearly correlate with the powerlaw flux, tracing a complete differ- 
ent pattern with the respect to type-A and -C QPOs. Type-B QPOs 
follow a different behavior also in a rms vs rms plane. In addition, 
they show a peculiar association to increases in count rate that could 
reflect changes in accretion rate and/or geometry in the system. 

A ll the types of QPOs can be explained through the precession 
model ( Ingram et a L 2009. [Ingram & Dondl2lnol. Ilngram & Pond 
l201lh as the result of the vertical Lense-Thirring precession of 
a optical translucent inner hot flow in a truncated disk geometry. 
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Table 3: Power-spectral classification and variability parameters. Only observations with evidence of low frequency QPOs are 
listed. 



# 


Obs. ID 


MJD 


Outburst 


Hardness ratio 


Total fractional rms 


QPO centroid Frequency 


QPO type 


State 


1 


701 09-0 l-07-00b 


52411.601 


2002 


0.228 ± 0.001 


7.9 ±0.1 


5.8 ±0.1 


B 


SIMS 


2 


70110-01-14-00 


52416.596 


2002 


0.252 ±0.001 


9.2 ±0.1 


6.4 ±0.1 


B 


SIMS 


3 


70110-01-15-00 


52419.238 


2002 


0.251 ±0.001 


9.0 ±0.1 


5.7 ±0.03 


B 


SIMS 


4 


70108-03-02-00b 


52419.432 


2002 


0.209 ± 0.001 


9.3 ±0.2 


5.6 ±0.1 


B 


SIMS 


5 


70110-01-47-00b 


52532.749 


2002 


0.222 ± 0.001 


7.2 ±0.2 


6.2 ±0.1 


B 


SIMS 


6 


70110-01-89-00 


52707.915 


2002 


0.256 ± 0.004 


9.7 ±0.7 


0.9 ±0.05 


B 


SIMS 


7 


90110-02-01-03 


53232.993 


2004 


0.257 ±0.001 


9.4 ±0.2 


4.1 ±0.04 


B 


SIMS 


8 


90704-01-02-00 


53233.389 


2004 


0.271 ±0.001 


9.1 ±0.2 


4.4 ±0.2 


B 


SIMS 


9 


60705-01-84-02 


53333.899 


2004 


0.254 ± 0.001 


8.8 ±0.1 


5.2 ±0.1 


B 


SIMS 


10 


91105-04-10-00 


53466.753 


2004 


0.293 ± 0.003 


12.8 ±0.7 


3.4 ±0.1 


B? 


SIMS 


11 


92035-01-04-00 


54147.011 


2007 


0.262 ± 0.001 


8.4 ±0.1 


6.7 ±0.2 


B 


SIMS 


12 


92085-01-03-01 


54162.665 


2007 


0.259 ±0.001 


8.6 ±0.1 


6.4 ±0.1 


B 


SIMS 


13 


92704-03-10-00 


54231.604 


2007 


0.264 ± 0.002 


6± 1 


1.0 ±0.1 


B 


SIMS 


14 


92704-03-10-11 


54232.658 


2007 


0.276 ± 0.002 


10.0 ±0.5 


1.7 ±0.1 


B 


SIMS 


15 


92704-03-10-12 


54233.565 


2007 


0.335 ± 0.002 


13.7 ±0.4 


1.8 ±0.2 


B? 


SIMS 


16 


95409-01-15-02 


55304.714 


2010 


0.256 ±0.001 


8.3 ±0.2 


5.6 ±0.1 


B 


SIMS 


17 


95409-01-15-06 


55308.983 


2010 


0.236 ±0.001 


6.7 ±0.1 


5.9 ±0.2 


B 


SIMS 


18 


95409-01-16-05 


55315.695 


2010 


0.262 ± 0.001 


9.3 ±0.1 


6.1 ±0.2 


B 


SIMS 


19 


95409-01-17-00 


55316.114 


2010 


0.262 ± 0.001 


9.3 ±0.1 


5.9 ±0.1 


B 


SIMS 


20 


95409-01-17-05 


55321.718 


2010 


0.243 ± 0.001 


7.5 ±0.1 


5.3 ±0.1 


B 


SIMS 


21 


95409-01-17-06 


55322.230 


2010 


0.238 ± 0.001 


7.5 ±0.2 


5.2 ±0.1 


B 


SIMS 


22 


95409-01-18-00 


55323.210 


2010 


0.259 ±0.001 


8.6 ±0.1 


5.5 ±0.1 


B 


SIMS 


23 


95335-01-01-07 


55324.189 


2010 


0.255 ±0.001 


8.7 ±0.2 


5.3 ±0.04 


B 


SIMS 


24 


95335-01-01-00 


55324.254 


2010 


0.247 ±0.001 


8.5 ±0.1 


5.3 ±0.03 


B 


SIMS 


25 


95335-01-01-01 


55324.393 


2010 


0.240 ± 0.001 


7.4 ±0.1 


5.1 ±0.02 


B 


SIMS 


26 


95335-01-01-05 


55326.175 


2010 


0.228 ±0.001 


7.1 ±0.2 


4.9 ±0.03 


B 


SIMS 


27 


95335-01-01-06 


55326.280 


2010 


0.223 ±0.001 


6.3 ±0.2 


4.9 ±0.03 


B 


SIMS 


28 


95409-01-18-04 


55327.041 


2010 


0.235 ± 0.002 


8.5 ±0.4 


4.8 ±0.1 


B 


SIMS 


29 


95409-01-18-05 


55327.262 


2010 


0.234 ± 0.001 


7.8 ±0.3 


4.9 ±0.04 


B 


SIMS 


30 


95409-01-19-00 


55330.300 


2010 


0.223 ±0.001 


6.5 ±0.6 


4.7 ±0.05 


B 


SIMS 


31 


r\s a r\r\ r\ i r\ a f\ A 

96409-01-04-04 


55585.947 


2010 


0.282 ± 0.003 


10.0 ± 0.6 


2.0 ± 0.1 


B 


SIMS 


32 


n^/inn m r\A nc 
y 04UV-U 1 -U4-U.) 


ooooD.oyo 


2010 


U.zby ± U.UU4 


y ± 2 


n n _l_ n i 
0.9 ± 0.1 


1 1 

D 




33 


96409-01-05-01 


55591.615 


2010 


0.292 ± 0.003 


10 ± l 


1.7 ± 0.1 


B 


SIMS 


34 


yo4uy-ui-u5-02 


OOO93.O02 


om n 
2010 


n odd _l n nno 
0.299 ± 0.003 


O ± 3 


i o _i_ n i 
1.8 ± 0.1 


D 


CTA /TO 


35 


701 09-0 l-07-00a 


52411.601 


2002 


0.213 ±0.001 


2.6 ±0.1 


7.0 ±0.5 


A 


SIMS 


36 


70108-03-02-00a 


52419.432 


2002 


0.214 ±0.001 


2.8 ±0.2 


6.7 ±0.5 


A 


SIMS 


37 


70108-03-02-00a2 


52419.432 


2002 


0.210 ±0.001 


2.8 ±0.2 


6.8 ±0.4 


A 


SIMS 


38 


70110-01-45-00 


52524.948 


2002 


0.222 ± 0.001 


2.2 ±0.2 


7.2 ±0.6 


A 


SIMS 


39 


70109-01-23-00 


52529.580 


2002 


0.209 ± 0.001 


2.9 ±0.2 


7.4 ± 1.3 


A 


SIMS 


40 


70109-01-24-00 


52536.358 


2002 


0.204 ± 0.001 


2.4 ±0.2 


8.0 ±0.9 


A 


SIMS 


41 


92085-01-02-06 


54160.896 


2007 


0.217 ±0.001 


2.6 ±0.2 


7.8 ±0.7 


A 


SIMS 


42 


92085-01-03-04 


54165.527 


2007 


0.210 ±0.001 


2.7 ±0.1 


7.7 ±0.7 


A 


SIMS 


43 


40031-03-02-05 


52388.054 


2002 


0.766 ± 0.003 


29.8 ±0.3 


0.20 ±0.01 


C 


LHS 


44 


70109-01-05-01G 


52391.318 


2002 


0.763 ± 0.003 


29.7 ±0.2 


0.22 ±0.02 


C 


HIMS 


45 


70109-01-06-00 


52400.83 


2002 


0.697 ± 0.002 


22.2 ±0.1 


1.26 ±0.01 


c 


HIMS 


46 


70108-03-01-00 


52400.853 


2002 


0.694 ± 0.002 


22.0 ±0.1 


1.30 ±0.01 


c 


HIMS 


47 


70110-01-10-00 


52402.492 


2002 


0.562 ± 0.002 


20.1 ±0.1 


4.20 ±0.08 


c 


HIMS 


48 


70109-04-01-00 


52405.58 


2002 


0.354 ±0.001 


15.45 ±0.04 


5.46 ±0.01 


c 


HIMS 


49 


70109-04-01-01 


52405.713 


2002 


0.356 ±0.001 


15.44 ±0.02 


5.45 ±0.01 


c 


HIMS 


50 


70109-04-01-02 


52406.07 


2002 


0.360 ± 0.001 


15.6 ±0.1 


5.34 ±0.02 


c 


HIMS 


51 


70110-01-11-00 


52406.701 


2002 


0.342 ± 0.001 


15.0 ±0.1 


5.82 ±0.02 


c 


HIMS 


52 


70110-01-12-00 


52410.528 


2002 


0.266 ± 0.001 


11.5 ±0.1 


8.1 ±0.2 


c 


HIMS 


53 


70109-01-37-00 


52694.922 


2002 


0.237 ± 0.002 


10.7 ±0.4 


8.6 ±0.2 


c 


HIMS 


54 


70128-02-02-00 


52696.355 


2002 


0.272 ± 0.001 


12.1 ±0.1 


8.02 ± 0.04 


c 


HIMS 


55 


50117-01-03-01 


52706.767 


2002 


0.364 ± 0.003 


19.2 ±0.3 


6.7 ±0.1 


c 


HIMS 


56 


50117-01-03-00 


52706.84 


2002 


0.363 ± 0.002 


18.7 ±0.2 


6.77 ± 0.02 


c 


HIMS 



Continued on next page 
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Table 3 - continued from previous page 



# 


Obs. ID 


MJD 


Outburst 


Hardness ratio 


Total fractional rms 


QPO centroid Frequency 


QPO type 


State 


57 


70109-02-01-00 


52709.859 


2002 


0.293 ± 0.002 


14.2 ±0.2 


8.0 ±0.1 


C 


HIMS 


58 


70109-02-01-01 


52709.991 


2002 


0.289 ± 0.002 


14.3 ±0.4 


8.1 ±0.1 


C 


HIMS 


59 


60705-01-56-00 


52710.715 


2002 


0.306 ± 0.003 


15.1 ±0.4 


7.8 ±0.1 


C 


HIMS 


60 


70110-01-94-00 


52724.226 


2002 


0.395 ± 0.004 


18.2 ±0.4 


6.1 ±0.1 


C 


HIMS 


61 


70110-01-95-00 


52727.252 


2002 


0.480 ± 0.005 


22.8 ±0.5 


4.7 ±0.1 


C 


HIMS 


62 


60705-01-59-00 


52731.562 


2002 


0.586 ± 0.003 


25.2 ±0.2 


2.9 ±0.1 


C 


HIMS 


63 


60705-01-68-00 


53218.11 


2004 


0.763 ± 0.005 


30.8 ±0.3 


0.5 ±0.0 


C 


LHS 


64 


60705-01-68-01 


53222.24 


2004 


0.728 ± 0.004 


26.9 ±0.2 


1.03 ±0.04 


C 


HIMS 


65 


60705-01-69-00 


53225.40 


2004 


0.714 ± 0.003 


24.9 ±0.2 


1.3 ±0.0 


C 


HIMS 


66 


90704-01-01-00 


53226.43 


2004 


0.705 ± 0.003 


24.9 ±0.1 


2.0 ±0.1 


C 


HIMS 


67 


60705-01-69-01 


53228.99 


2004 


0.646 ± 0.003 


23.7 ±0.2 


2.9 ±0.2 


C 


HIMS 


68 


60705-01-70-00 


53230.96 


2004 


0.437 ± 0.002 


18.3 ±0.2 


4.3 ±0.1 


C 


HIMS 


69 


90110-02-01-02 


53232.34 


2004 


0.386 ± 0.002 


17.5 ±0.1 


5.2 ±0.1 


C 


HIMS 


70 


90110-02-01-00 


53232.40 


2004 


0.361 ±0.001 


16.35 ±0.04 


5.8 ±0.1 


C 


HIMS 


71 


90704-01-11-00 


53472.33 


2004 


0.584 ± 0.004 


24.9 ±0.3 


2.7 ±0.1 


C 


HIMS 


72 


92035-01-02-01 


54133.922 


2007 


0.771 ± 0.003 


30.8 ±0.2 


0.28 ±0.01 


C 


LHS 


73 


92035-01-02-02 


54135.033 


2007 


0.771 ± 0.003 


30.8 ±0.2 


0.30 ±0.01 


C 


LHS 


74 


92035-01-02-03 


54136.015 


2007 


0.766 ± 0.003 


30.1 ±0.2 


0.37 ±0.01 


C 


LHS 


75 


92035-01-02-04 


54136.997 


2007 


0.759 ± 0.003 


29.8 ±0.2 


0.43 ±0.01 


C 


LHS 


76 


92035-01-02-08 


54137.851 


2007 


0.748 ± 0.003 


27.5 ±0.3 


0.55 ±0.02 


C 


HIMS 


77 


92035-01-02-07 


54138.83 


2007 


0.731 ± 0.002 


25.8 ±0.2 


0.90 ±0.01 


C 


HIMS 


78 


92035-01-02-06 


54139.942 


2007 


0.686 ± 0.002 


22.4 ±0.1 


0.99 ±0.01 


C 


HIMS 


79 


92035-01-03-00 


54140.204 


2007 


0.670 ± 0.002 


21.7 ±0.1 


1.13 ±0.01 


C 


HIMS 


80 


92035-01-03-01 


54141.055 


2007 


0.621 ± 0.002 


20.7 ±0.1 


1.68 ±0.01 


C 


HIMS 


81 


92035-01-03-02 


54142.036 


2007 


0.547 ± 0.002 


19.6 ±0.1 


2.45 ±0.01 


C 


HIMS 


82 


92035-01-03-03 


54143.019 


2007 


0.461 ± 0.002 


18.3 ±0.1 


3.52 ±0.01 


C 


HIMS 


83 


92428-01-04-00 


54143.870 


2007 


0.411 ±0.001 


17.1 ±0.1 


4.34 ± 0.02 


c 


HIMS 


84 


92428-01-04-01 


54143.951 


2007 


0.419 ±0.001 


17.2 ±0.1 


4.23 ±0.02 


c 


HIMS 


85 


92428-01-04-02 


54144.086 


2007 


0.424 ± 0.002 


17.4 ±0.1 


4.13 ±0.03 


c 


HIMS 


86 


92428-01-04-03 


54144.871 


2007 


0.380 ± 0.001 


16.4 ±0.1 


4.99 ± 0.03 


c 


HIMS 


87 


92035-01-03-05 


54145.114 


2007 


0.343 ± 0.001 


14.9 ±0.1 


5.80 ± 0.03 


c 


HIMS 


88 


92085-01-03-00 


54161.669 


2007 


0.295 ± 0.001 


12.9 ±0.1 


7.1 ±0.1 


c 


HIMS 


89 


92085-01-03-02 


54163.698 


2007 


0.288 ± 0.001 


12.5 ±0.1 


7.3 ±0.2 


c 


HIMS 


90 


92085-01-03-03 


54164.557 


2007 


0.296 ± 0.001 


12.4 ±0.1 


7.0 ±0.2 


c 


HIMS 


91 


92704-03-11-00 


54234.839 


2007 


0.493 ± 0.004 


21.3 ±0.4 


4.0 ±0.1 


c 


HIMS 


92 


92704-03-11-01 


54235.791 


2007 


0.549 ± 0.005 


22.1 ±0.4 


3.3 ±0.1 


c 


HIMS 


93 


92704-04-01-01 


54236.446 


2007 


0.577 ± 0.004 


22.8 ±0.6 


3.0 ±0.2 


c 


HIMS 


94 


92704-04-01-02 


54236.513 


2007 


0.587 ± 0.004 


24.8 ±0.4 


2.7 ±0.1 


c 


HIMS 


95 


92704-03-12-00 


54236.591 


2007 


0.598 ± 0.005 


25.5 ± 1.1 


2.7 ±0.4 


c 


HIMS 


96 


92704-04-01-04 


54237.356 


2007 


0.593 ± 0.004 


25.2 ±0.5 


2.8 ±0.1 


c 


HIMS 


97 


92704-04-01-05 


54237.421 


2007 


0.584 ± 0.004 


24.5 ±0.4 


3.1 ±0.3 


c 


HIMS 


98 


92704-03-12-01 


54237.488 


2007 


0.596 ± 0.004 


24.8 ±0.3 


2.7 ±0.1 


c 


HIMS 


99 


95409-01-12-04 


55286.727 


2010 


0.787 ± 0.003 


32.4 ±0.4 


0.22 ±0.01 


c 


LHS 


100 


95409-01-13-03 


55288.367 


2010 


0.783 ± 0.003 


31.5 ±0.3 


0.2 ±0.1 


c 


LHS 


101 


95409-01-13-00 


55289.618 


2010 


0.777 ± 0.003 


31.0 ±0.3 


0.26 ±0.01 


c 


LHS 


102 


95409-01-13-04 


55290.722 


2010 


0.781 ± 0.003 


31.7 ±0.2 


0.29 ±0.01 


c 


LHS 


103 


95409-01-13-02 


55291.649 


2010 


0.775 ± 0.003 


31.6 ±0.3 


0.32 ±0.01 


c 


LHS 


104 


95409-01-13-05 


55292.779 


2010 


0.777 ± 0.003 


30.9 ±0.4 


0.38 ±0.02 


c 


LHS 


105 


95409-01-13-01 


55293.088 


2010 


0.772 ± 0.003 


30.6 ±0.4 


0.38 ±0.05 


c 


LHS 


106 


95409-01-13-06 


55294.124 


2010 


0.770 ± 0.003 


30.1 ±0.4 


0.47 ±0.02 


c 


LHS 


107 


95409-01-14-01 


55296.248 


2010 


0.734 ± 0.003 


26.7 ±0.3 


1.04 ±0.01 


c 


HIMS 


108 


95409-01-14-02 


55297.87 


2010 


0.672 ± 0.002 


21.9 ±0.1 


1.25 ±0.01 


c 


HIMS 


109 


95409-01-14-03 


55298.70 


2010 


0.648 ± 0.003 


21.6 ±0.2 


1.59 ±0.01 


c 


HIMS 


110 


95409-01-14-06 


55299.766 


2010 


0.564 ± 0.002 


20.0 ±0.1 


2.43 ±0.01 


c 


HIMS 


111 


95409-01-14-04 


55300.336 


2010 


0.563 ± 0.002 


20.0 ±0.2 


2.38 ±0.01 


c 


HIMS 


112 


95409-01-14-07 


55300.923 


2010 


0.515 ±0.002 


19.4 ±0.1 


2.92 ±0.0 


c 


HIMS 


113 


95409-01-14-05 


55301.789 


2010 


0.454 ± 0.002 


17.7 ±0.2 


3.64 ±0.02 


c 


HIMS 


114 


95409-01-15-00 


55302.196 


2010 


0.425 ± 0.002 


17.5 ±0.2 


4.15 ±0.03 


c 


HIMS 


115 


95409-01-15-01 


55303.604 


2010 


0.346 ± 0.001 


14.6 ±0.1 


5.65 ± 0.04 


c 


HIMS 
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Table 3 - continued from previous page 



# 


Obs. ID 


MJD 


Outburst 


Hardness ratio 


Total fractional rms 


QPO centroid Frequency 


QPO type 


State 


116 


95409-01-17-02 


55318.441 


2010 


0.308 ± 0.001 


13.7 ±0.1 


6.67 ±0.19 


C 


HIMS 


117 


96409-01-06-01 


55598.700 


2010 


0.491 ± 0.004 


26.7 ± 1.5 


4.52 ±0.30 


C 


HIMS 
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Table 4: Columns are: observation number, reduced x 2 , inner disc temperature (kT), inner disc radius R (assuming a distance 
of 10 kpc and an inclination of 30 °), photon index T, fold Energy E/ ;<2 (corresponding to high energy cutoff), total flux, hard 
flux and disk flux calculated in the 2 - 20 keV band and expressed in units of Crab flux. 



# 


reduced \ 


T " Innradius 


p 1 nn 


r 


E/oZd 


Ftot/FCrab 


^hard^FCrab 


^diskfFCrab 


instruments 


1 


0.92 


octroi 


43t 3 


2.6±g;i 


- 


0.501 


0.063 


0.321 


PCA+HEXTE 


2 


1.01 


o.93^n 4 4 


41±S 


2.6l°? 


- 


0.462 


0.070 


0.270 


PCA+HEXTE 


3 


1.11 


0-93±°-gS 


38±S 


9 «+ - 2 
^• 8 -0 2 


- 


0.434 


0.066 


0.235 


PCA+HEXTE 


4 


0.81 


0-90^3 


43l« 


2.61°'? 


- 


0.395 


0.048 


0.259 


PCA+HEXTE 


5 


0.96 


0-911^3 


45l 4 


3-0±g"l 


- 


0.494 


0.060 


0.297 


PCA+HEXTE 


6 


0.81 


o.6itg:S 


33+? 1 


o O+0.2 
d -°-0.3 


- 


0.027 


0.004 


0.016 


PCA+HEXTE 


7 


0.85 


o.8ot°;°§ 


37l 4 


9 7+0.2 
z - ' -0.4 


- 


0.178 


0.027 


0.104 


PCA+HEXTE 


8 


0.84 


0.85t°:£ 
0.90+°;°| 


31+ 3 


2.61°,:! 


- 


0.175 


0.029 


0.099 


PCA+HEXTE 


9 


1.07 


38l 3 


2-5l°-l 


- 


0.314 


0.047 


0.195 


PCA+HEXTE 


10 


0.93 


O-81-o 07 


19±S 


2-4l°1 


- 


0.056 


0.011 


0.030 


PCA+HEXTE 


11 


0.92 


oq+0 04 

u.yo_Q Q2 


46l 3 


9 7+0.I 
— 03 


- 


0.653 


0.099 


0.367 


PCA+HEXTE 


12 


1.08 


Q2+ 03 
o-y^_o 03 


44l 4 


9 7+0.1 
^* ' —0 1 


- 


0.536 


0.080 


0.307 


PCA+HEXTE 


13 


0.89 


o.76t°;°S 


25l 4 


9 9+O.3 
^•^—0.3 


- 


0.054 


0.008 


0.035 


PCA+HEXTE 


14 


0.99 


0.781°"°* 


22i 5 3 


2 l+°- 3 

^■J-— 2 


- 


0.053 


0.009 


0.033 


PCA+HEXTE 


15 


1.04 


0-8±g-l 


16+ 3 


2 l+°' 3 


- 


0.046 


0.010 


0.025 


PCA+HEXTE 


16 


1.26 


0-96l°° 3 4 


40+ 4 


9 ,1+0.04 
^•^_0 04 


- 


0.496 


0.071 


0.318 


PCA+HEXTE * 


17 


1.18 


o.9olg-S 


48+ 4 


9 r+0.1 
^• J — 05 


- 


0.492 


0.062 


0.329 


PCA+HEXTE * 


18 


1.00 


n q 4 +o-03 

O.»1_o 04 


38+ 4 


9 r + 0.1 
z ' d — 1 


- 


0.438 


0.068 


0.261 


PCA+HEXTE * 


19 


1.25 


o.9ii°:S 


43l| 


9 o+O.l 
z ' J -0.04 


- 


0.432 


0.064 


0.275 


PCA+HEXTE * 


20 


1.23 


0.93l°;° 4 


39l| 


2 4+ ' 1 


- 


0.382 


0.052 


0.252 


PCA+HEXTE * 


21 


1.54 


o.89l° ;° 4 


44+ 4 


9 oe+0.04 
z - ou — 0.04 


- 


0.383 


0.049 


0.263 


PCA+HEXTE * 


22 


1.74 


0-93l°;° 4 


38l 4 


9 qo+0.04 
z - oo -0 .04 


- 


0.384 


0.057 


0.242 


PCA+HEXTE * 


23 


1.07 


0.89l°;° 3 3 


41±t 


9 7+0.4 
z - ' -0.4 


- 


0.391 


0.056 


0.225 


PCA+HEXTE * 


24 


1.07 


0.87l°°t 


43l 4 


9 7+0.2 
^•'—0 2 


- 


0.385 


0.054 


0.226 


PCA+HEXTE * 


25 


1.12 


0-86±g;°* 


45l 3 


9 7+0.2 
z - ' —0.2 


- 


0.380 


0.051 


0.231 


PCA+HEXTE * 


26 


1.03 


0.84«-° 4 3 


47 +5 


2 9+°' 2 

^• y _0 4 


- 


0.371 


0.046 


0.219 


PCA+HEXTE * 


27 


1.28 


0.85l° ° 4 


46 + s 


9 7+O.3 
^•'—0 3 


- 


0.357 


0.043 


0.230 


PCA+HEXTE * 


28 


1.26 


0.89±g"gS 


4ll 5 5 


9 o+O.l 
1 


- 


0.327 


0.042 


0.230 


PCA+HEXTE * 


29 


1.50 


0.901°$ 


4011 


2 4+ - 1 

^•^_0 03 


- 


0.342 


0.045 


0.232 


PCA+HEXTE * 


30 


1.21 


0.86l° ° 4 


4611 


24 +o.i 

Zl - 4: — 1 


- 


0.330 


0.039 


0.232 


PCA+HEXTE * 


31 


1.27 


n 71 +0.06 
u ' ' ± -0 06 


3ll^ 4 


2 +0 - 2 


- 


0.055 


0.009 


0.038 


PCA+HEXTE * 


32 


1.78 


94+ - 14 

O-!"_0.03 


13±1 


i-4±g;i 


- 


0.039 


0.006 


0.029 


PCA+HEXTE * 


33 


0.94 


+0.06 
u - ' ' -0.07 


22l 4 


i-8±g;i 


- 


0.042 


0.007 


0.028 


PCA+HEXTE * 


34 


1.31 


0.85l°;° 6 4 


1511 


i-8±g;i 


- 


0.035 


0.007 


0.023 


PCA+HEXTE * 


35 


1.09 


0-94l°°4 


44l« 


2 6+ ' 1 


- 


0.488 


0.055 


0.336 


PCA+HEXTE 


36 


1.31 


0.90±g"g? 


44 +3 3 


2 6+ ' 1 


- 


0.394 


0.045 


0.271 


PCA+HEXTE 


37 


1.31 


0.901°,°,* 


44l 3 


2 6+ ' 1 


- 


0.394 


0.045 


0.271 


PCA+HEXTE 


38 


0.97 


0-95l°g 4 4 


441^ 


o 7+0.I 
^•'—0 2 


- 


0.531 


0.065 


0.342 


PCA+HEXTE 


39 


0.85 


0.951^3 


43l 3 


2.6i°-; 


- 


0.480 


0.053 


0.335 


PCA+HEXTE 


40 


0.89 


0.90l°, : ° 3 


48l 3 


2.6l° : i 


- 


0.448 


0.047 


0.311 


PCA+HEXTE 


41 


0.86 


0.90l°° 4 


51±2 


9 7+0.I 
z - ' -0.1 


- 


0.535 


0.060 


0.357 


PCA+HEXTE 


42 


1.08 


0.88+°-° 4 


49l 3 


9 S+ ' 1 
2-8-0.2 


- 


0.461 


0.050 


0.306 


PCA+HEXTE 


43 


1.19 




1 74+0.02 

'^-0.02 


119+f 


0.401 


0.236 


0.000 


PCA+HEXTE 


44 


0.88 






i 70+O.02 
1- '8-0.02 


133+1° 


0.407 


0.238 


0.000 


PCA+HEXTE 


45 


0.91 






1 Q6+ ' 02 


122+ 11 


0.430 


0.230 


0.000 


PCA+HEXTE 


46 


1.50 






1 qo+0.02 
i - y °-0.02 


1291° 


0.428 


0.227 


0.000 


PCA+HEXTE 


47 


0.84 






9 4l -+0.03 
z -*°-0.03 




0.410 


0.173 


0.000 


PCA+HEXTE 


48 


0.80 


1 oi +0 - 05 

i - ul -o.o5 


26+ 3 


2 59+0- 07 

^.O»_ 22 




0.494 


0.122 


0.164 


PCA+HEXTE 


49 


1.12 


n 99+0-04 

u - aa -0.04 


27l 3 


9 r fi +0.06 
Z - OD -0.15 




0.499 


0.123 


0.171 


PCA+HEXTE 


50 


0.85 


1 0+ 01 


24l 3 


9 7+0.1 
z - ' -0.1 




0.471 


0.119 


0.140 


PCA+HEXTE 


51 


0.84 


1 qo+0.06 
1 - uo -0.06 


28+ 3 


9 9+O.4 
^•^-0.2 




0.460 


0.106 


0.215 


PCA+HEXTE 


52 


0.94 


94+ - 04 


39l 4 


9 r + 0.1 
z - d -0.2 




0.467 


0.078 


0.259 


PCA+HEXTE 


53 


0.89 


69+ 05 


30±1 


9 4 +o.6 

^•^-0.2 




0.046 


0.006 


0.030 


PCA+HEXTE 


54 


1.36 


r, 7 O+0.03 

u - ' °-0.03 


23+| 


Z - D -0.1 




0.041 


0.007 


0.022 


PCA+HEXTE 


55 


0.67 


0.9i° ;l 




9 r+0.2 




0.034 


0.009 


0.011 


PCA+HEXTE 
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Table 4 - continued from previous page 



# 


reduced x 2 




*DUk 


r 




Ftot/Fcrab 


Fhard/Fcrab 




instruments 


56 


1.11 




15+* 


2.5^ 


- 


0.035 


0.009 


0.013 


PCA+HEXTE 


57 


0.97 




19±1 


9 9 +0.2 
—0 2 


- 


0.032 


0.006 


0.019 


PCA+HEXTE 


58 


0.74 


0.7±g : l 


27±l 


2A+° i 


- 


0.035 


0.006 


0.021 


PCA+HEXTE 


59 


1.01 


0-71;;!! 


211? 


2-8l°1 


- 


0.026 


0.005 


0.013 


PCA+HEXTE 


60 


1.29 






o 9+0.1 
°- — 1 


- 


0.035 


0.010 


0.000 


PCA+HEXTE 


61 


1.03 




- 


2 9+ ' 1 


- 


0.034 


0.012 


0.000 


PCA+HEXTE 


62 


1.17 


_ 


_ 


2 48+ ' 03 

—0 06 


_ 


0.033 


0.015 


0.000 


PCA+HEXTE 


63 


1.02 




- 


1 oq+0.03 
J--OO_ 03 


- 


0.102 


0.061 


0.000 


PCA+HEXTE 


64 


1.15 




- 


l-98t°o : °o| 


- 


0.126 


0.070 


0.000 


PCA+HEXTE 


65 


1.20 




- 


1-97±S"S? 


- 


0.132 


0.073 


0.000 


PCA+HEXTE 


66 


1.40 




- 


2 01 + ' 03 

o 02 


- 


0.136 


0.074 


0.000 


PCA+HEXTE 


67 


0.90 




- 


9 r>o+0.03 

^•^0_ Q Q3 


- 


0.134 


0.066 


0.000 


PCA+HEXTE 


68 


0.94 


1 2+ ' 1 

— 1 


8+? 


o q+0.1 
° — 1 


- 


0.124 


0.040 


0.038 


PCA+HEXTE 


69 


1.09 


i 17 +0.09 

— 09 


11+? 


1 7+ ' 7 

^* ' — 1 


6411? 


0.127 


0.035 


0.065 


PCA+HEXTE 


70 


0.88 


i -,0+0.05 
J---LO_ 0.08 


1211 


y — 2 


IO2I38 


0.127 


0.033 


0.064 


PCA+HEXTE 


71 


1.14 




- 


2 4+ ' 1 

—0 1 




0.046 


0.021 


0.000 


PCA+HEXTE 


72 


1.30 


- 


- 


i 09+0.02 

i-O^-O 02 


1631H 


0.367 


0.214 


0.000 


PCA+HEXTE 


73 


1.19 


- 


- 


1.80lH 2 2 


139+J 1 


0.371 


0.216 


0.000 


PCA+HEXTE 


74 


1.52 






oq+0.02 

±.O0_Q Q2 


14611 3 


0.386 


0.222 


0.000 


PCA+HEXTE 


75 


1.25 


_ 


_ 


1.84l° : ° 2 2 


137+* 1 


0.413 


0.236 


0.000 


PCA+HEXTE 


76 


1.17 


- 


- 


1-86+To 2 . 


13lll« 


0.424 


0.239 


0.000 


PCA+HEXTE 


77 


1.13 


- 


- 


1.92l°'° 2 2 


140+" 


0.433 


0.238 


0.000 


PCA+HEXTE 


78 


1.03 


- 


- 


9 Qcr + 0.02 
z - UJ -0.02 


142+1? 


0.434 


0.223 


0.000 


PCA+HEXTE 


79 


1.05 






2 09+ ' 02 

i>uy -0 02 




0.436 


0.219 


0.000 


PCA+HEXTE 


80 


0.81 


- 


- 


9 9 r+0.02 
z - zd -0.02 




0.443 


0.206 


0.000 


PCA+HEXTE 


81 


0.92 






9 rqn+0-004 
z-.uoy_Q 004 




0.475 


0.191 


0.000 


PCA+HEXTE 


82 


1.00 


1 +tMJ 


19 +6 3 


2.46±° : ^ 


- 


0.501 


0.166 


0.085 


PCA+HEXTE 


83 
84 


1.00 
0.96 


0.980+°'go4 

±.UOO_q ggtj 


23At°i 

is.8ir 2 


9 coQ+0.003 

Z-.rJOf _q QQ3 

9 ^07+0.004 

Z.UO * _Q QQ4 


- 
- 


0.533 
0.529 


0.154 
0.157 


0.121 
0.103 


PCA+HEXTE 
PCA+HEXTE 


85 


0.85 


1 02+ 00 


i9.o±8i 


9 fi oo+0.004 

^.UOO_Q QQ4 


- 


0.522 


0.157 


0.098 


PCA+HEXTE 


86 


0.88 


1.0241^04 


22.6l°i 


2.69910-^4 


- 


0.538 


0.141 


0.142 


PCA+HEXTE 


87 


1.10 


0.98l° °02 


30±i 


2.631?"?? 


_ 


0.554 


0.127 


0.203 


PCA+HEXTE 


88 
89 


1.05 
1.26 


n 929+ 000 

KJ.UZii? _Q QQ2 

922 +o: °° 


on r> + 0.2 
OO.Z q 2 

38.5lf 2 


9 roo+0.003 
z.joz_q 003 
9 r 9 r + 0:003 
a.ozjO_q 003 


_ 
- 


0.488 
0.462 


0.092 
0.085 


0.246 
0.240 


PCA+HEXTE 
PCA+HEXTE 


90 
91 


0.97 
0.98 


0.898^q'oo2 


39.61^ 


2 641 +0 - 003 

u — 003 
9 OO+0.02 


- 
- 


0.465 
0.048 


0.088 
0.018 


0.222 
0.000 


PCA+HEXTE 
PCA+HEXTE 


92 


1.32 


- 


- 


9 rq+0.02 
^■00_ 01 


- 


0.046 


0.019 


0.000 


PCA+HEXTE 


93 


1.23 


- 


- 


2.46+Toi 


- 


0.045 


0.020 


0.000 


PCA+HEXTE 


94 


1.28 


- 


- 


2.45l? : ?l 


- 


0.045 


0.021 


0.000 


PCA+HEXTE 


95 


1.30 


- 


- 


9 qO + 0.02 


- 


0.045 


0.021 


0.000 


PCA+HEXTE 


96 


1.16 


- 


- 


2-45l -°oi 


- 


0.045 


0.021 


0.000 


PCA+HEXTE 


97 


1.25 


- 


- 


9 r 9 + 0.01 


- 


0.044 


0.020 


0.000 


PCA+HEXTE 


98 


1.97 


- 


- 


9 45+O.O6 




0.043 


0.020 


0.000 


PCA+HEXTE 


99 


1.47 


_ 


_ 


-, 01+0.02 

i-Oi_0 02 


_ 


0.285 


0.168 


0.000 


PCA+HEXTE * 


100 


1.77 


- 


- 


1 8n+ ' 02 


- 


0.293 


0.172 


0.000 


PCA+HEXTE * 


101 


1.55 






1 09+0.02 

i-O^-O 03 




0.300 


0.175 


0.000 


PCA+HEXTE * 


102 


1.86 


- 


- 


1 09+0:03 


- 


0.307 


0.179 


0.000 


PCA+HEXTE * 


103 


1.73 






1 or+0.02 
J-SO-O.OS 




0.317 


0.183 


0.000 


PCA+HEXTE * 


104 


1.68 






1 oc+0.02 




0.317 


0.183 


0.000 


PCA+HEXTE * 


105 


1.24 






1 Q1+0.03 
1 - 81 -0.03 


171.893l 2 g 


0.335 


0.194 


0.000 


PCA+HEXTE * 


106 


1.36 






1 or+0.03 
1 - oo -0.02 




0.340 


0.194 


0.000 


PCA+HEXTE * 


107 


1.08 






1 - az -0.04 


191.624! 3 ! 


0.347 


0.190 


0.000 


PCA+HEXTE * 


108 


1.25 






2 10+ 003 

z - lu -0.03 


0.367 


0.184 


0.000 


PCA+HEXTE * 


109 


1.50 






9 9n +0.00 
^•^ u -0.00 




0.375 


0.179 


0.000 


PCA+HEXTE * 


110 
111 


1.17 
1.10 


n 7+0.57 
u - ' -0.3 

o.8i° ;| 


22+1 88 


9 2 q+0.09 
^•^ a -0.16 
9 9+O.2 
^•^-0.2 




0.405 
0.392 


0.160 
0.158 


0.043 
0.046 


PCA+HEXTE * 
PCA+HEXTE * 


112 
113 


1.55 
1.19 


0.981 ,; 22 

1 n97 +0.004 
l.u^l _o.004 


20t£° 
20.9±g;l 


9 1+0.2 

9 -I or + 0.004 
Z.lOO_ 004 


153.151^° 


0.392 
0.405 


0.144 
0.129 


0.083 
0.122 


PCA+HEXTE * 
PCA+HEXTE * 
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Table 4 - continued from previous page 



# 


reduced x 1 




*DUk 


r 




Ftot/Fcrab 


Fhard/Fcrab 




instruments 


114 
115 
116 
117 


1.19 
1.40 
1.11 
1.24 


0.9i- 

1.0021°;™ 

n qko+0.002 
u.ooo_ 002 


30+f 
28.9+° \ 
42.8±g;l 


9 97/1+0.004 
^•^'1_0.004 
9 /IK7+0-004 
Z -^J' -0.004 

2.95316l°,;J° 




0.429 
0.454 
0.387 
0.032 


0.132 
0.100 
0.074 
0.012 


0.137 
0.214 
0.197 
0.000 


PCA+HEXTE * 
PCA+HEXTE * 
PCA+HEXTE * 
PCA+HEXTE * 
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